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Abstract 

We develop a renormalized continuum field theory for a directed poly- 
mer interacting with a random medium and a single extended defect. The 
renormalization group is based on the operator algebra of the pinning poten- 
tial; it has novel features due to the breakdown of hyperscaling in a random 
system. There is a second-order transition between a localized and a delocal- 
ized phase of the polymer; we obtain analytic results on its critical pinning 
strength and scaling exponents. Our results are directly related to spatially 
inhomogeneous Kardar-Parisi-Zhang surface growth. 
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Low-dimensional manifolds in media with quenched disorder are objects en- 
countered in a large variety of different physical realizations. Obvious examples 
are interfaces in disordered bulk media and random field systems |], 3 or magnetic 
flux lines in dirty superconductors 0, but there is also a deep connection to the 
problem of nonequilibrium surface growth j|, [j| and randomly driven hydrodynam- 
ics || . Furthermore, the theory serves as a simple paradigm for more complicated, 
fully frustrated random systems such as spin glasses 0. 

A one dimensional manifold in a random medium is a phenomenological con- 
tinuum model for a (single) magnetic flux line in type-II superconductors with 
impurities |3], [T3J, where the flux line interacts with an ensemble of quenched point 
defects (represented by a random potential). In addition to these point impuri- 
ties there may also be extended (e.g columnar or planar) defects in the system. 
Experiments that systematically probe the effect of this kind of impurities have 
recently become possible in high temperature superconductors ||. The statistics 
of the line configurations is governed by an energetic competition: point defects 
tend to roughen the flux line; it performs large transversal excursions in order to 
take advantage of locally favourable regions. An attractive extended defect, on the 
other hand, suppresses these excursions and, if it is sufficiently strong, localizes the 
line to within a finite transversal distance The two regimes are separated by 
a second order phase transition where the localization length £j_ diverges. In con- 
trast to temperature-driven transitions, it involves the competition of two different 
configuration energies rather than energy and entropy and is hence governed by a 
zero-temperature renormalization group fixed point. 

The system is described by an effective Hamiltonian 

n= i dt -v^+po^y w 

Here r(t) denotes the displacement vector of the flux line (also called directed poly- 
mer) in d' transversal dimensions, as a function of the longitudinal "timelike" co- 
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ordinate t. The random potential V(r, t), Gauss-distributed with V(r,t) = and 



V(r, t)V(r', t') = 2cr 2 S d '(r — r')5(t — £'), models the quenched point disorder. Av- 
erages over disorder are denoted by an overbar, thermal averages by brackets (...). 
The last term in the Hamiltonian describes the interaction with the extended de- 
fect. In this letter, we concentrate on columnar defects, $(£) = 8 d '(r(t)), but most 
of the results can straightforwardly be generalized to planar defects. 

This model gains additional interest since it is related to nonequilibrium critical 
phenomena |S]]. If Z(r,t) denotes the restricted partition sum over all paths ending 
at a fixed given point (r,t), the "height field" h(r,t) = log Z(r, t) obeys the 
evolution equation 

§ = v^h + ^{Vhf +V- Po 5 d \r) (2) 

with uq = (2/3) -1 and Ao = 1. This is the Kardar-Parisi- Zhang (KPZ) equation of 
directed surface growth [IIJ with an additional term describing a local inhomogene- 
ity in the rate of mass deposition onto the surface. 

Quite a few authors have studied these models. Early numerical work for d' = 1 



indicates a derealization transition at a finite defect strength Pq \ H[ |T2"f . In other 
large scale simulations [ 13| , however, it is found that arbitrarily weak defects local- 
ize the polymer for d' = 1, but a finite defect strength is necessary for d' > 1. This 
result is supported by an approximate renormalization treatment for the problem 
on a hierarchical lattice fl3|], by scaling arguments |L4|, and by an approximate 



functional renormalization |T5[ . In |L6[ a Wilson type renormalization is discussed, 
but its consistency is unclear. All of these approaches rely on nonsystematic approx- 
imations; and the status of the transition has remained controversial. The problem 
has so far defied attempts at an exact solution even for d' = 1, in contrast to the 
related problem of disorder- induced depinning from a rigid wall |TT|, |T7 . 



This letter is devoted to a field-theoretic study of the derealization transition. 
The zero-temperature fixed point of directed lines in a random potential (with 
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Po = 0) is a scale-invariant continuum field theory. Its two basic exponents, the 
roughness exponent ( and the anomalous dimension of the disorder- averaged free 
energy —u, are defined in eqns. @ and (|4]) below. This theory is the starting 
point for a systematic perturbation theory in the pinning potential; it involves 
an e-expansion with borderline dimension d' — 1. In contrast to standard cases 
like 4 -theory, here even the unperturbed system is a field theory with complicated 
multipoint correlation functions, due to the non-thermal averaging over the disorder. 
Nevertheless, two fundamental properties of the local pinning field can be 
obtained in terms of the exponents ( and u: (i) its scaling dimension and (ii) 
the form of its operator product expansion, see eqns. (||) and (||) below. These 
properties determine the renormalization group equations for the pinning strength 
to leading order, and hence the phase diagram of the system. We find a transition 
at po — for d' < 1, and at a finite (nonuniversal) pinning strength Pq for d' — 2. 
The renormalization may equivalently be carried out in the framework of the KPZ 
dynamics, eq. This is discussed at the end of this letter, together with some 
implications for inhomogeneous growth processes. 

To derive the renormalization group equations for the depinning problem, it is 
necessary to study the polymer configurations in the unperturbed random system 
(i.e. for po = 0). The leading scaling behavior on large scales of disorder-averaged 
correlation functions is due to the sample-to-sample fluctuations of the polymer 



paths of minimal energy ||18|| . Typical transversal excursions of the paths, given e.g. 
by the two-point function 

A 2 (t! - t 2 ) = ((r^-r^)) 2 ) ~ \h - t 2 \ 2 < , (3) 

are characterized by the roughness exponent (. It is larger than for thermal fluctu- 
ations, namely ( = 2/3 for d! = 1 and ( ~ 5/8 for d' = 2 (see e.g. and references 
therein). The exponent — u is the anomalous dimension of the disorder-averaged 
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free energy F = —f3 1 logTr exp(— (37i), whose universal part has the scaling form 

F(T, R) ~ T w F{R/T C ) (4) 
in a finite system of transversal size R and longitudinal size T. Due to a "tilt 



symmetry", these exponents obey the relation u = 2£ — 1 (see e.g. [18|, [19], |5|). 

This asymptotic scaling is determined by a renormalization group fixed point 
where the temperature (3~ l is an irrelevant scaling variable, which nevertheless 
turns out to be crucial to what follows. It governs the crossover from the Gaussian 
fixed point (a 2 = po — 0), that describes the free thermal scaling on scales much 
smaller than the crossover length £y , to the disorder fixed point. The finite-size 
free energy (|]) has the crossover scaling form F(T, (3~ l ) = T w jF(T/£n) at fixed 
R/T^. In the infrared, F is independent of in the ultraviolet, F ~ /3 _1 T°. 
Comparing these two limits gives ~ ^ , i.e. —to is the anomalous dimension 
of the temperature at the disorder fixed point. In an analogous way, the two-point 
function (|3]) has the scaling form A 2 (t,/3 -1 ) = t 2 ^ J(f3~ l t~^) with a temperature- 
independent infrared limit and the ultraviolet asymptotic behavior A 2 (t, (3~ l ) = 
/? _1 /f: given by the Gaussian theory. Hence the corrections to scaling of ([3]) due to a 
small temperature (which are illustrated in fig. 1) are given by 

A^/T 1 ) = t 2 ^ + cf3- 1 t 2 ^ + ... (5) 

with 2£ — u — 1, assuming analyticity of the scaling function J . 

Now we discuss the correlation functions of the local pinning field at the 



disorder fixed point. The one-point function ($(£)) gives the probability that at 
time t, the polymer is at the origin r = 0, averaged over thermal and disorder 
fluctuations. In a system of infinite longitudinal size, but finite transversal size R 
(and periodic boundary conditions), one has by translational invariance 



($(t)) = Lr x with x = d'(, (6) 
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where L is a longitudinal scale defined by L = R 1 '^. In the renormalization dis- 
cussed below, L has the role of an infrared cutoff, and it generates the renormal- 
ization group flow. The exponent x is the scaling dimension of the field $ at the 
disorder fixed point. 



The full multipoint correlation functions (3>(ti) • • -$(t m )) give the probability 
that a (single) path crosses the line r = at given times ti . . .t m . To discuss 
the short distance properties of these objects, specifically consider the two-point 



function (<&(ti)Q(ti)) for \t\— t 2 | «C L. In this limit, it depends on the infrared cutoff 
as L~ x , i.e. in the same way as the one-point-function (Q). Hence asymptotically, 



it factorizes into ($(ti)) and the L-independent "return probability" to the origin 
(which is simply the inverse spread of the paths ~ A(t)~ d ' given by (0)), 

WhWM ~ 1*1 - t 2 \- d ' c (i - cp-% -t 2 r+.. (7) 

Again the leading singularity is due to the sample-to-sample fluctuations of the 
minimal energy paths, while the correction term is due to the thermal fluctuations 
around these paths. The leading, temperature-independent singularity equally oc- 



curs in the thermally disconnected two-point function ($(£1)) ($(t 2 )) |2IJ. Hence 



in the connected two-point function $(t2)) c , only the subleading singularity 

survives. An analogous argument applies to the singularities in any correlation 



function (. . . $(£)$(£') . . .) as \t — t'\ 0. Therefore the relation 

$(t)$(t / )~ c /r 1 |t-t / r x ~ w $(t) (8) 

(with a constant c > 0) is valid as an operator identity, i.e. inserted into an arbi- 



trary connected correlation function (. . . $(£)$(£') . . .) c . The notion of an operator 
algebra that encodes the universal short- distance properties of correlation functions 
is familiar in field theory |ET|. The new feature of (^) is that the leading singularity 



is governed by a correction-to-scaling exponent. This is a consequence of the break- 
down of hyperscaling at the zero-temperature disorder fixed point and has direct 
implications for the renormalization of the pinning problem to which we now turn. 



The renormalized perturbation theory for the Hamiltonian ([I]) is constructed 
along the lines of ref . f22[ . The universal part of the disorder-averaged free energy 
density / = lim-r^oo-F 1 /TR d ' can be expanded in powers of the dimensionful bare 
coupling constant po, 
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f(po, L) - /(0, L) = p- l R- d ' ' T Jm , (9) 



where 

J m = J dt 2 ■ - • dt m ($(0)$(t 2 )---$(t m )) c . (10) 

A weak defect potential distorts the minimal energy paths of the unperturbed sys- 
tem, the dominant paths reorganize exploiting the low-lying excitations. As dis- 
cussed above, the statistics of these excitations is encoded in the connected corre- 
lation functions that appear in J m . The leading ultraviolet singularities of these 
integrals are dictated by the operator algebra (§). In an analytic continuation to ar- 
bitrary d', these show up as poles in e(d') = 1 — (d'( + uj), which serves as expansion 
parameter. Inserting ([]) and (g) into (^), we find to second order 

Jip^T) - WI) = -L"- 1 ($) Lq - ^y^o 2 ) + O (w 3 , e°w 2 ) ■ (11) 

where wq = p§U is the dimensionless bare coupling constant. The pole in e(d') 
can be absorbed into a renormalized coupling constant w = Z(w)wq with Z(w) = 
1 — (c/e)w + O {w 2 ). Its renormalization group flow 

Ld L w = e{d')w ~cw 2 + (V) (12) 

determines the large scale behavior of the perturbed system. For e(d') > 0, i.e. for 
d! < 1, the perturbation is relevant: for any attractive bare defect potential, the 
renormalized coupling is driven towards large attractive values. Hence the flux line 
is localized by an arbitrary weak attractive columnar defect. The localization length 
diverges as £j_ ~ \po\ ~ u± with u± = ((d')/e(d') when the defect strength approaches 
zero from below. In the borderline dimension d' — 1 an attractive defect potential is 



marginally relevant: the line is still localized by an arbitrary weak columnar defect, 
but with an essential singularity in the localization length £j_ ~ exp(2/3c|w|). 

For s(d!) < 0, i.e. for dl > 1, a weak defect is an irrelevant perturbation. The 
transition to a pinned state now takes place at a finite critical strength pQ (which 
however depends on the microscopic scales of the system and is hence nonuniversal). 



It is governed by the nontrivial fixed point w* = e(d')/c < of (|12"D . Close to 
the transition, the localization length diverges as £j_ ~ |po — Po\ u± y where u±_ = 
C(d')/y*(d') and y*(d') is given by the ^-expansion y* = —e(d') + 0(e 2 ). 

Additional insight into this problem may be gained by the mapping onto the 
KPZ equation (0). From this stochastic equation, one constructs in a standard 
way the generating functional Tr exp(— S[h, h]) of the dynamical correlation func- 
tions (denoted by ((. . .))) in terms of the height field h and the "conjugate" field h 
[p5fl . Insertions of this field generate response functions, e.g. ((h(r,t) Uih(ri,ti))) = 
{{5 h(r, t)/Hi5 V(ri, U))). In the dynamical action, there is a term Si = po J dth(0, t), 
the analogue of the pinning term in ([!]). The perturbation series (§) for the re- 



stricted free energy (3 1 logZ(r, t) = {{h(r,t))) is in one-to-one correspondence with 
the dynamical perturbation series; we have 



pm-i ($(ti) . . . m m) y = ^ Q h ( Q) h y.. h ( Q> tm))) (13) 

Hence h(0, t) = /3$(t) is a field of scaling dimension x = d' ( + u; by virtue of (H), 
it obeys the short-distance algebra 

h(0, t)h{0, t') ~ c \t - t'\~* h{0, t') + ... . (14) 

Its leading singularity is no longer a correction-to-scaling exponent; the peculiarity 
of the correlation functions written in terms of the h fields is rather that they have to 
be computed in the nontrivial "vacuum" state h(r,t). This makes the nonunitarity 
of this theory manifest, which is generated by the averaging over disorder. 

For po < 0, the term Si describes an excess mass deposition onto the growing 
surface at r = 0. This term breaks the translational invariance. For stationary 



growth ((h(r,t))) = v t + H(r) in a one- dimensional system of size R, it results in 
an excess growth velocity and an approximately triangular surface profile H(r) f9j. 
From the mapping to the polymer system, one concludes that for £±(w) <C R, the 
excess velocity scales as (v(w, R) - v(0,R)) ~ ~ exp(2/3cw). The same 

essential singularity shows up in the slope \dH/dr\ ~ (v(w,R) — v(0, R)) 1 / 2 . The 
response function has the form {{h(r,t)h(ri,ti))) = for t — > oo and 

r, £j_ -R. All of these quantities are accessible in numerical simulations which 
could provide a useful test of the results discussed in this letter. 

In summary, we have shown that a class of field theories with quenched random- 
ness shares with conventional field theories the notion of a short- distance algebra 
of its scaling operators, whence they are amenable to systematic renormalization. 

Note added: In the final stages of this work we received a preprint by Hwa and 
Nattermann on the depinning problem. They exploit the mapping onto the 
KPZ equation (0). In d! = 1, they obtain the algebra fll4|) by using a mode- 
coupling approximation for the dynamical correlation functions. 

We thank J. Krug, R. Lipowsky, and L.-H. Tang for useful discussions. 
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Figure 1 

Ensemble of paths with a common starting point in a typical disordered sample. 
The thick line denotes the path of minimal energy; its transversal fluctuations 
are given by A 2 (t) ~ t 2 ^ . The thin lines represent thermal fluctuations of width 
S 2 (t) = (r 2 (t)) c ~ j3~H around the minimal path. 
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